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1 Introduction: relativistic harmonic-oscillator problem 

The simplest and perhaps most straightforward generalization of the Schrodinger operators of 
standard nonrelativistic quantum theory towards the inclusion of relativistic kinematics leads 
to Hamiltonians H that involve the relativistic kinetic energy, or relativistically covariant form 
of the free energy, of a particle of mass m and momentum p, given by the square-root operator 



and a coordinate-dependent static interaction potential V(x). In the one-body case, they read 



The eigenvalue equation of this Hamiltonian is usually called the "spinless Salpeter equation." 
It may be regarded as a well-defined approximation to the Bethe-Salpeter formalism [1] for the 
description of bound states within relativistic quantum field theories, obtained when assuming 
that all bound-state constituents interact instantaneously and propagate like free particles [2]. 
Among others, it yields semirelativistic descriptions of hadrons as bound states of quarks [3,4]. 

In general, the above semirelativistic Hamiltonian H is, unfortunately, a nonlocal operator: 
either the relativistic kinetic energy, T(p), in configuration space or, in general, the interaction 
potential in momentum space is a nonlocal operator. Because of the nonlocality it is somewhat 
difficult to obtain rigorous analytical statements about the solutions of its eigenvalue equation. 
Thus sophisticated methods have been developed to extract information about these solutions; 
for details and comparisons of the various approaches, consult, for example, the reviews [5-10]. 

Analytical or, at least, semianalytical 1 expressions for both upper and lower bounds on the 
eigenvalues of some self-adjoint operator may be found by combining the minimum-maximum 
principle [11-13] and appropriate operator inequalities [14-18]. The outcome of this procedure 
is sometimes called the "spectral comparison theorem." In Sec. 3 below, we will use this kind of 
bounds in order to estimate the accuracy of our findings for the eigenvalues of the operator (1). 
Accordingly, we recall in Appendix A the proof of the "translation" of some inequality satisfied 
by two operators into the corresponding relations between their discrete eigenvalues, by briefly 
sketching all basic assumptions and the line of argument. A very systematic path for obtaining 
such operator inequalities is provided by rather simple geometrical considerations summarized 
under the notion "envelope theory" [19-24]. These envelope techniques may be generalized to 
systems composed of arbitrary numbers of relativistically moving interacting particles [25-27]. 
For particular potentials V, semianalytical lower bounds on the ground-state energy eigenvalue 
of the semirelativistic Hamiltonian (1), and hence on the entire spectrum of H, can be found by 
the appropriate generalization of the local-energy theorem [13,28-30] to our case of relativistic 
kinematics [23,31], or by applying the optimized (Beckner-Brascamp-Lieb) version of Young's 
inequality for convolutions to some integral formulation of the spinless Salpeter equation [32]. 

Purely numerical solutions of the spinless Salpeter equation may be computed in numerous 
ways. The semirelativistic Hamiltonian H can be approximated by some effective Hamiltonian 
which is of nonrelativistic shape but uses parameters that depend on expectation values of the 
momenta [33,34]. Upper bounds of, in principle, arbitrarily high precision on the eigenvalues of 
a self-adjoint operator can be found [18,35-39] with the help of the Rayleigh-Ritz (variational) 

1 We regard a bound on some eigenvalue of a given self-adjoint operator as semianalytical if it can be derived 
by the — in general, numerical — optimization of an analytically known expression over a single real variable. 





(1) 



2 



technique as immediate consequence of the minimum-maximum theorem [11-13]. The spinless 
Salpeter equation may also be converted into an equivalent matrix eigenvalue problem [40-45]. 
A particular role for H is played by the spherically symmetric harmonic-oscillator potential 

V(x) = ar 2 , r = \x\ , a > ; 

this potential defines the "relativistic harmonic-oscillator problem," posed by the Hamiltonian 

H= ^p 2 + m 2 + ar 2 . (2) 

The eigenvalue equation of H, for eigenstates H\ip) = E involves only one parameter: 
upon factorizing off some overall energy scale a 1 / 3 by performing the canonical transformation 

and rescaling both mass m and eigenvalue E according to m — a 1 / 3 /j, and E = a 1//3 e, it reads 

( vV + ^ + r 2 ) M=eM • 
In momentum-space representation this eigenvalue equation reduces to a Schrodinger problem 

[-A p + V(p)]il>(p) = eil>(p) , (3) 
with an interaction potential reminiscent of the square root of the relativistic kinetic energy: 

v( P ) = vW/' 2 • 

We would like to take advantage of this facet of the relativistic harmonic-oscillator problem 
in order to derive for its bound-state eigenfunctions, in this analysis, a compact expression and 
to move thereby beyond only numerically calculated exact solutions without having to rely on 
perturbation theory, maybe paving the way for the eventual construction of analytic solutions. 

2 Analytical solutions for bound-state eigenfunctions 

We exploit the fact that, in contrast to the general case, for a harmonic-oscillator potential (as, 
with due care, for any potential of the form V oc r 2n , n = 1, 2, 3, . . .) the eigenvalue equation of 
the Hamiltonian H is an ordinary differential equation, parametrized by // and its eigenvalue e. 

Focusing on the ground state or purely radial excitations, let us introduce, for eigenstates 
with vanishing relative orbital angular momentum £, the reduced radial wave function y(p) by 

y{p) = V^p^{p) (^ = 0), 
which is, of course, subject to the normalization condition 

oo 

J dp\y(p)\ 2 = 1 • 



In accordance with Eq. (3) the reduced radial wave function satisfies a reduced radial equation: 



%p) = [V( P )-e}y(p) 



(4) 
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From its definition, this reduced radial wave function y(p) has to vanish at the origin: y(0) = 0. 
Moreover, the analysis of the normalizable solutions of the eigenvalue equation (4) reveals that 
y(p) behaves like p for small p, that is, for p< 1; hence, its derivative with respect to p at the 
point p = is a nonvanishing constant, which may be absorbed into the overall normalization: 



1 . 



x ^ P" 



We construct all solutions of Eq. (4) in form of Taylor-series expansions by using the ansatz 

yip) 

with the expansion coefficients 

d n ?/ 

0,1,2,... 



C n = -T— (0) , 



&p r - 



n 



The solution of the eigenvalue equation (4) is then clearly tantamount to the determination of 
the expansion coefficients c n . The first three of these expansion coefficients are known trivially: 



c = 2/(0) = , 

* - >-» 

d 2 y 



c 2 = 



dp 



2 (0) = [(V-e)y](0) = 0. 



[For the sake of notational simplicity, we suppress, in accordance with our above remark, in the 
following that normalization factor of y(p) which guarantees its unity norm and assume y(p) to 
be normalized such that the value of the first nonvanishing expansion coefficient is one: c\ — 1.] 
Upon insertion of the eigenvalue equation (4) followed by the application of Leibniz's theorem, 
the nontrivial expansion coefficients c n ,n > 3, may be shown to satisfy the recurrence relation: 



c n+2 



dp n+2 



(0) 



fc=0 



dp n 
d k {V 



dp 



n 



= (fi-e)c n + J2 
k=i 

n 

= {fi-e)c n + Y, 



k=i 



d 2 y 


<o) = ! 


dp 2 




-e) 


d n ~ k y 


k 


dp n ~ k 


I) 




I) 





d n [(V-e)y] 



dp r - 



(0) 



(0) 



n = 1,2,3, 



(5) 



Here, for the last step, we abbreviated the k-th derivative of the potential V by a coefficient d k : 
d k V t 



d k = H 



dp 1 



■(0) 



d k Vx 2 + 1 
dx k 



(0), 



_ P 
x = — 



H> , k = 0,1,2,... . 
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In the case /x — m — 0, the solutions involve Airy's function Ai(z); cf., e. g., Refs. [18,23,33,34]. 
According to the above definition, we have do — 1. Furthermore, by inspection of the function 
f(x) = \/x 2 + 1, it is easy to convince oneself that all odd derivatives of f(x) vanish at x — 0, 



d 2fc+1 = for all k = 0, 1, 2, . . . , (6) 
whereas all even derivatives of /(x) at x = necessarily satisfy the (simple) recurrence relation 

d 2k+2 = (1-4 k 2 )d 2k , A; = 0,1,2,... . 
By induction, the solution of this recurrence relation for the nonvanishing coefficients d 2 k reads 



d 



(2fc-l) 



(2k-2)\ 



1 2 



fc = 1,2,3,... . 



(7) 



_2 fe - 1 (A;-l)! 

Taking into account the observation (6), we obtain C3 = fj,—e and, for the coefficients c n , n > 4, 
c n +2 = (n ~ e) c n + y 2 k J ^ ^ d2h ° n ~ 2k ' n = 2, 3, 4, . . . , 



where 



n 
2 

n — 1 



for n even, n = 2, 4, 6, . . . 



for n odd, n — 3, 5, 7, . . . 



Thus the recurrence relation (5) for all expansion coefficients c n decomposes into one involving 
only the even coefficients c 2n , n = 2, 3, 4, ... , and one involving only the odd coefficients c 2n+ i, 
n — 2, 3, 4, . . . ; recalling Co = and c 2 = 0, we conclude that all the even coefficients c n vanish: 

c 2n = for all n — 0, 1, 2, . . . . 

With the result (7), the recurrence relation for the (nonvanishing) odd coefficients finally reads 



c 2n +3 = - e) c 2n +i + 



k=l 



(/i - e) c 2n+1 + 

k=i \ 

n = 1,2,3, ... . 



2n + 1 
2k 

2n + 1 
2A; 



,l-2fe 



<^2fc C2n- 



2fc+l 



^fc+i/i 1 " 2 ' {-\f-\2k-l) 



k-i 



(2k-2)\ 
2 k - 1 (k-l)\ 



(8) 



In summary, upon constructing the relevant expansion coefficients according to this recurrence 
relation the analytical expressions for all reduced radial wave functions (of £ = bound states) 
y{p) of the relativistic harmonic-oscillator problem (2) are given by the power-series expansion 

00 p2n+l 



y(P)= E C ™ (2» + l)|- 



(9) 



The various solutions y(p) of the eigenvalue equation (3) are characterized or discriminated 
by different sets of expansion coefficients c n . By construction, apart from the first coefficient c\, 
all expansion coefficients for a given solution depend on the corresponding energy eigenvalue e. 
Evaluating the recurrence relation (8) for just the first term, the series (9) explicitly starts with 



yip) = p 1 1 + (a* - e) |j + 
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3 Explicit bound-state eigenfunctions and energy levels 

The central result of our present investigation of the "relativistic harmonic-oscillator problem" 
defined by the Hamiltonian H of Eq. (2) is an analytical expression for the reduced radial parts 
y(p) of all eigenfunctions ip(p) of H for vanishing angular momentum in form of a Taylor series: 

y{p) = VIirpip(p) = C 2n+1 777- - TT, , 

71=0 \ Z L )- 

where the corresponding expansion coefficients C2„+i, n — 0,1,2, ... , are either fixed to c\ — 1, 
C3 = fi—e or, for C2„+i, n > 2, determined by the recurrence relation (8). For a given solution of 
this eigenvalue problem these expansion coefficients and thus the resulting reduced radial wave 
function y(p) depend on the parameter /x = m/a 1 ^ 3 and the corresponding energy eigenvalue e. 
Taking into account the necessary requirement of normalizability of Hilbert-space eigenstates, 
the inversion of the latter relation may be exploited to determine the energy eigenvalues e from 
the knowledge of the dependence of e on the coefficients c n , as derived in the preceding section. 
In order to fulfil its normalization condition, any solution y(p) must vanish in the limit p — > 00: 

limy(p) = 0. (10) 

Fixing the energy eigenvalue e of a chosen bound state in this manner and using this particular 
value of the parameter e in the expansion (9) then yields the corresponding wave function y{p). 

Our principal concern is beyond doubt the semianalytical approach developed in Sec. 2 and 
summarized above. Nevertheless, it might be instructive to construct explicitly a few examples 
of solutions in numerical or graphical form. These results can be compared with the outcome of 
some straightforward (but merely numerical) integration of the Schrodinger equation (4). This 
will provide a useful check of the correctness of our solutions and justify the present formalism. 

In actual computations, the infinite series (9) has to be truncated, for practical purposes, to 
a reasonably large but definitely finite number N of terms considered in this expansion of y{p): 

N p 2n+l 

yW "£ C2 " +1 (2^TT)!' 

In this case, the wave function y{p) will approach zero, as required by the constraint (10), not 
at infinity but already for a finite value, say p, of the momentum p — before it starts to diverge. 
This technique gives the energy eigenvalues e with a precision determined by one's choice of N: 
e = e(N). Likewise, the momentum boundary p will also change with N: p = p(N). For a given 
value of \i, which quantifies the relative importance of particle mass m and harmonic-oscillator 
coupling strength a, every wave function resulting from this truncation procedure involves two 
dimensionless parameters: the relevant energy eigenvalue e(N) of the Hamiltonian (2), and the 
characteristic momentum p(N). Their values will be determined simultaneously, in accordance 
with the above requirement on y(p) [to approach zero at p(N)] by an appropriate fit procedure. 
In other words, the value of p(N), in particular, cannot be varied freely; it is fixed for chosen N. 

Let us illustrate this procedure for both ground state and first radial excitation, i. e., for the 
two bound states defined by vanishing orbital angular momentum and radial quantum number 
n T = 0, 1, resp. Figure 1 shows the corresponding reduced radial wave functions y{p) for /i — 30 
as obtained by inspecting the functional form of y(p) resulting from different choices of e and p, 
if taking into account 45 terms in their Taylor series (9), that is, if choosing iV = 44 in Eq. (11). 
Table 1 summarizes the relevant numerical parameter values emerging from such construction. 
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1.75 r 




(b) 

Figure 1: Radial eigenfunctions in momentum space y(p) = \^Att pi/j(p) of the ground state (a) 
and the first radial excitation (b) of the relativistic harmonic-oscillator problem, defined by the 
Hamiltonian H = ^p 2 + /i 2 +r 2 , with /i = 30, and for 45 terms in the Taylor expansion of y(p). 

The exact results can be easily computed with the aid of a (standard) integration technique 
designed for solving the Schrodinger equation numerically [46] . For our two examples the exact 
wave functions y{p) prove to be practically indistinguishable, at least by the eye, from the ones 
extracted from the series (11) with N = 44. This explains why we refrain from plotting also the 
former in Fig. 1. For the momentum range depicted in Fig. 1, that is, for < p < p, the relative 
differences of the areas under corresponding curves are (of the order) 10~ 8 ; more precisely, they 
are given by 2x 1CT 8 for n r = 0, the ground state, and 3x 1CT 8 for n r = 1, the first excited level. 
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Table 1: Dimensionless (by scaling) energy eigenvalues e(N) and characteristic momenta p(N) 
of the relativistic harmonic-oscillator problem posed by the Hamiltonian H = y/p' 2 + /x 2 +r 2 as 
well as the classical turning points p of the corresponding nonrelativistic motion for the ground 
state and the first radial excitation (identified by their radial quantum number n r = 0, 1), with 
mass- vs. -coupling strength parameter pi = 30, and N = 44 or N = 14 in our Taylor series (11). 



N 


n r 


e(N) - pi 


p(N) 


p 


44 





0.38627 


10 


4.82 




1 


0.89864 


12 


7.36 


14 





0.38854 


8.5 


4.82 



From a straightforward consideration of the "classical turning points" of the corresponding 
(and well understood) nonrelativistic harmonic-oscillator problem defined by the Hamiltonian 

#nr = /i + — + r 2 , 
2 pi 

the maximum "classical" momenta, p, are found, in terms of the radial quantum number n r , as 

p 2 = (4n r + 3)y / 2/I , n r = 0, 1, 2, . . . . 

By inspecting Table 1 we note with satisfaction that for both energy levels under consideration 
the numerical values of the suitable p turn out to be far beyond their classical counterparts p. 

In order to get, at least, some vague idea of the dependence of our findings on the amount of 
truncation represented by N < oo, we inspect the ground-state wave function y{p) constructed 
again for pi = 30 but by truncating the expansion (9) to the rather modest number of 15 terms, 
which means to set N = 14 in Eq. (11). Figure 2 confronts this approximate wave function y(p) 
with its exact behaviour for the ground state. While for small and intermediate momenta there 
is still perfect agreement with the exact result [46] , we observe a clearly discernible discrepancy 
between approximate and exact curve for large momenta. Table 1 tells us that even for N = 14 
our crucial momentum p is still comfortably above the corresponding classical turning point, p. 
Moreover, comparing the cases N = 14 and iV = 44, we learn that the value of p increases with 
increasing number N. Of course, the naive expectation would be that p behaves like p — > oo for 
N — > oo, that is, when removing the truncation and restoring the full series expansion for y{p). 

The minimum number of terms to be taken into account in the Taylor-series expansions (9) 
required in order to achieve some given precision of one's results will depend, of course, on both 
the bound state under study and the desired accuracy. From our above remarks we feel entitled 
to conclude that a reasonable (and manageable) number N ~ 40 produces satisfactory results. 

To our knowledge, at present the best semianalytical upper and lower bounds to the energy 
eigenvalues of the relativistic harmonic-oscillator problem are provided simultaneously by the 
combination of minimum-maximum principle with operator inequalities [18] and the envelope 
theory [19,21]: at least for the relativistic harmonic oscillator the envelope bounds [19,21] may 
be shown [23] to be quantitatively equivalent to the bounds derived in Appendix A of Ref. [18]. 
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Figure 2: Momentum-space wave function y(p) of Fig. 1(a) resulting from consideration of only 
15 terms in its polynomial approximation (11) [full line], in comparison with the corresponding 
exact ground-state wave function of the relativistic harmonic-oscillator operator [dashed line] . 



However, a discussion, in full generality, of all implications of such operator inequalities for the 
eigenvalues of the operators considered appears clearly off the mainstream of our presentation. 
Therefore, as already promised in the Introduction, the general relationship is demonstrated in 
Appendix A. The bounds we need here are derived from this general theorem by specializing to 
the case of the relativistic harmonic-oscillator problem posed by the Hamiltonian operator (2); 
all operator inequalities required by this procedure may be generated by, e. g., envelope theory. 
For a bound state of vanishing relative orbital angular momentum £, that is, for a purely radial 
excitation, identified by the radial quantum number n T — 0, 1, 2, . . . (identical to the number of 
nodes of the corresponding wave function), the bounds on the dimensionless eigenvalue e read 




min I \l a 2 H ^ + r 2 I < : < uiin 

r>0 \ V r 2 / r>0 

where, in three spatial dimensions, our envelope-theory upper-bound parameter Pu is given by 

3 

Pu = 2n r + - , n, = 0,1,2, ... , 

while the lower-bound parameter Pl required for our envelope bounds is related to the zeros zq 
of Airy's function Ai (z) [47] (-z = 2.338107, 4.087949, 5.520560, 6.786708, 7.944134, . . .) by 

Pl = 2 



(^) , Ai(, ) = 0. 



Resulting values of Pl for the lowest-lying £ = bound states are listed in Table 2 [9,19,21-23]. 
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Table 2: Numerical values of the parameter Pl used by envelope theory for the lower bounds on 
the energy levels of the relativistic harmonic-oscillator problem in three spatial dimensions, for 
the lowest-lying £ = bound states identified by their radial quantum number n r = 0, 1, 2, . . . . 




1 

2 
3 
4 



1.3760835 
3.1813129 
4.9925543 
6.8051369 
8.6182269 



Table 3 compares, for the lowest four £ = energy levels (identified by their radial quantum 
number n T — 0, 1, 2, 3) of our Hamiltonian (2), the approximate eigenvalues, s(N), obtained by 
the present approach by a truncation of the power series (11) to N — 14 or N = 44 terms, with 
the corresponding semianalytical upper (en) and lower (ei) bounds mentioned above as well as 
with the (numerically exact) eigenvalues e num , computed by a method developed for the purely 
numerical solution of (nonrelativistic) Schrodinger equations [46] . For N = 14, the polynomial 
in e resulting from the suitably adapted boundary condition (10) has only two real roots at all. 
Moreover, both of these approximate values are still above our (semianalytical) upper bounds. 
In contrast to this rather crude approximation, for N = 44 the eigenvalues e(N) of already the 
three lowest energy levels fit perfectly to the ranges spanned by the semianalytical bounds. For 
the ground state, in particular, e(44) reproduces the exact result at least to five decimal places. 

Table 3: "Compact-origin" eigenvalues e(N), their upper (en) and lower (e L ) bounds, and their 
exact values e num for the lowest £ = states of the Hamiltonian H/a 1 ^ in Eq. (2), with /i = 30. 



Radial Excitation n r 







1 



e L - n 

^num (A 

e(N = 14) 
e(N = 44) 



H 
/i 



0.38668 
0.35478 
0.38627 
0.38854 
0.38627 



0.90032 
0.81862 
0.89857 
0.93616 
0.89864 



1.41179 
1.28222 
1.40768 



1.92111 
1.74440 
1.91366 



1.41032 1.94319 



4 Summary and Conclusions 

Our compact result for the reduced £ = eigenfunctions of the Hamiltonian (2) is given by the 
power series (9), with expansion coefficients c\ — 1, C3 = fJ-—£, and C2„+i, n > 2, determined by 
the recurrence relation (8). Both the numerical determination of all energy eigenvalues and the 
explicit construction of the corresponding eigenfunctions of the relativistic harmonic-oscillator 
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problem is then achieved by forcing our solutions to satisfy, in addition, the constraint imposed 
by the requirement of normalizability of bound-state wave functions. Comparing these explicit 
solutions with the outcomes of purely numerical integration procedures reveals that at least for 
the lowest-lying energy levels our semianalytical approach reproduces, already for a truncation 
of the Taylor series (9) to a moderate number of expansion terms, the exact solutions with high 
accuracy. Of course, if one is interested only in numerical solutions of the problem under study, 
their straightforward computation with the help of some integration algorithm should produce 
the desired result more easily than their extraction from our Taylor series by means of Eq. (10). 
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A Combination of minimum— maximum principle with 
operator inequality: "spectral comparison theorem" 

It is a simple exercise to relate discrete eigenvalues of two operators satisfying some inequality. 

Consider for some operator A, with domain T>(A), its eigenvalue equation A \a k ) = a k \aik), 
k — 0, 1, 2, . . . , for its set of eigenstates {\a k ), k = 0, 1, 2, . . .}, corresponding to its eigenvalues 



and likewise for some operator B, with domain T>(B), its eigenvalue equation B \/3 k ) = b k \f3 k ), 
k — 0, 1, 2, . . . , for its set of eigenstates {\j3 k ), k — 0, 1, 2, . . .}, corresponding to its eigenvalues 



Assume that both these operators A and B are self-adjoint: A^ = A, B^ = B. This implies that 
all their eigenvalues are real: a* k = a k , b* k — b k , k — 0, 1, 2, . . . . Let these eigenvalues be ordered 
according to a < ai < a 2 < ■ ■ • , b < b\ < b 2 < ■ ■ ■ . Consider only the discrete eigenvalues a k 
of A below the onset of the essential spectrum of the operator A. Assume that the operator A is 
bounded from below. Assume that the operators A and B are related by an operator inequality 
of the form A < B, which implies that B too is bounded from below. In order to derive, for any 
k — 0, 1, 2, . . . , the relationship between a k and b k , focus on some arbitrary (A;+l)-dimensional 
subspace D^+i of the domain T>(A) of A: D k+ i C T>(A). Employing the appropriate form of the 
minimum-maximum principle, the operator inequality A < B translates into an upper bound 
on the eigenvalue a k of A which involves all expectation values of B within this subspace Du+i- 



a k = 



(a k \ A \a k ) 
(a k \a k ) 



k = 0,1,2,... , 



h = 



(Pk\B\f3 k ) 



k = 0,1,2,... . 



a k < sup 

l^)eD fc+1 



(4>\A\4>) 



< sup 



(4>\B\ij) 



for all k = 0,1,2,... . 



(12) 




Now, in order to relate the supremum over the expectation values of B to the eigenvalues b k 
of B, consider a particular subspace D k+1 , namely, that space that is spanned by the first k+1 
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eigenvectors of the operator B, that is, by precisely those eigenvectors of B that correspond to 
the first k + 1 eigenvalues b , h, . . . , bk of B: D k+1 C V(B) C V(A). Then, clearly, every |V) in 
is a linear combination of the eigenstates i = 0, 1, . . . , k} of B, with coefficients cf. 

k 

W = Ec i |/3 i ) for all |V>) £ -Dfc+i • 

i=0 

For any subspace -Dfc+i, k — 0,1,2, ... , use of this expansion of |V) yields for its norm squared 

k 

(VIV) = E M 2 (AIA) for all |V) G £ fc+ i 

i=0 

and, with this and 6, < b k for alH = 0, 1, . . . , k, an upper bound on all expectation values of B: 

k k 

(V| B |V) = £ |q| 2 6, (AIA) <^E|q| 2 (AIA) = (VIV) for all |V) e , 

i=0 i=0 

which means 

(VI -B I 



(VIV) 

01 5 IV) 



< 6 fe for all |V) e Am , 
= 6 fc for |V) = |/3 fc ) G A+i . 



(VIV) 

Therefore the supremum of all expectation values of B over D k+1 is just the eigenvalue b k of B: 

sup ^77^ = 6 fc for all k = 0, 1, 2, . . . . 

|^>GD fc+1 (VIV) 

Thus, inserting this identity in the chain of inequalities (12) proves that corresponding discrete 
eigenvalues a k , b k of semibounded self-adjoint operators A, B that satisfy A < B are related by 

a k < b k for all k — 0, 1, 2, . . . . 
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